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IMTRQDUCTION 

The  purpose  of  this  paper  Is  to  examine  and  correlate  the 
various  physical  and  chemical  relations  proposed  to  explain  the  phenom- 
ena of  detonation.  In  particular,  the  detonation  of  solid  explosives. 
In  order  to  assure  a  reasonably  comprehensive  treatise,  It  Is  desirable 
at  the  outset  to  devote  some  attention  to  the  correction  of  certain 
widespread  misconceptions  regarding  explosives  and  their  behavior.  Two 
distinct  types  of  explosion  processes,  deflagration  and  detonation,  are 
evidenced  by  materials  referred  to  as  explosives.  It  has  been  common 
practice,  especially  In  the  early  literature  on  the  subject,  to  distin- 
guish between  the  two  on  the  basis  of  chemical  composition  of  the 
explosive  and  Its  Intended  application.  Deflagration  Is  Invariably 
associated  with  "low"  explosives  or  propellants  and  the  continual  employ- 
ment of  certain  explosives  such  as  nitrocellulose  and  black  powder  for 
their  propulsive  effect  has  led  to  the  erroneous  association  of  chemical 
compositions  with  the  type  of  behavior.  Similarly  the  detonation  process 
is  associated  with  primary  and  "high"  explosives  and  likewise  is  con- 
sidered to  be  the  sole  type  of  behavior  evidenced  by  other  chemical 
substances.  Unfortunately  the  above  classification  is  not  borne  out  by 
facts.  On  the  contrary,  experimental  results  indicate  that  all  explosives 
are  capable  of  exhibiting  both  modes  of  behavior.  In  view  of  this  it  is 
necessary  to  establish  scane  criterion  other  than  chemical  composition  to 
enable  one  to  distinguish  clearly  between  the  two  processes. 

The  deflagration  process,  in  its  broadest  sense,  can  be  con- 
sidered a  rapid  burning  in  which  the  essential  rate-determining  factor 
is  the  speed  of  the  chemical  reaction.  That  is,  the  rate  at  which  the 
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deflagration  proceeds  is  related  directly  to  the  kinetics  of  the  reaction 
taking  place  in  the  burning  zone. 

The  propagation  of  deflagration  is  closely  associated  vdth 
transport  mechanisms;  either  the  transport  of  heat  or  of  atoms  and  free 
radicals  acting  as  chain-carriers  or  a  combination  of  the  two.  As  long 
as  mechanisms  such  as  these  control  the  propagation,  there  is  established 
a  limiting  velocity  beyond  which  the  rate  of  deflagraticsi  cannot  pass. 
Lxtensive  investigations  have  been  carried  out  in  connection  with  the 
process  of  deflagration  in  a  gaseous  medium  and  several  excellent  treat- 
ments are  available  (1>2). 

In  the  case  of  detonation,  the  characteristic  featxire  distin- 
guishing the  process  is  the  appearance  of  a  pressure  wave  propagating  at 
a  high  velocity.  The  rate  at  which  detonation  pzM>ceeds  is  calculable  by 
means  of  hydrodynamic  and  thermodynamic  equations  without  the  necessity 
of  reaction-kinetic  data.  The  application  of  the  classical  hydrodynamic- 
thermodynamic  theory  requires  no  assumptions  regarding  chemical  mechanisms 
other  than  the  postulate  that  the  reactions  occur  at  such  a  rate  so  as 
not  to  be  a  limiting  factor. 

This  theory,  which  is  reviewed  in  detail,  is  based  upon 
ejqpressions  for  the  continuity  of  mass,  momentum  and  energy  across  a 
so-cailled  detonation  front  and  upon  the  assumption  that  the  stable 
detonation  velocity  is  the  minimum  one  compatible  with  the  conservation 
equations.  When  supplemented  by  an  equation  of  state  and  thermal  data 
for  the  products  of  reaction,  the  theory  permits  the  calculation  of 
detonation  velocity,  temperature  and  pressure,  as  well  as  other  parameters 
of  interest. 
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INITIATIQM  AND  PROPAGATION  OF  DETONATION 

The  currently  accepted  hydrodynamic  detonation  theory,  which 
>dll  be  discussed  in  a  later  section,  unfortunately  provides  no  informa- 
tion concerning  the  mechanism  by  vdiich  a  stable  detonation  wave  arises* 
It  is  well  established  by  experience  that  detonation  can  be  brought  about 
by  a  wide  variety  of  conditions,  including  mechanical  shock,  high 
temperature,  friction,  and  adiabatic  compression  of  interstitial  gas* 
Only  recently,  however,  have  successful  attempts  been  made  towards  the 
establishment  of  what  might  be  considered  a  generalized  theory  of  ini- 
tiation. The  work  of  S.  J.  Jacobs  and  D.  P*  HacDougall  (3)  at  the 
Explosives  Research  Laboratory  of  the  N.D.R.C.  involving  the  controlled 
deflagration  of  readily  detonable  substances  and  that  of  F*  P.  Bowden 
and  his  associates  (4)  on  the  isolation  of  factors  responsible  for 
thenoal  ignition  provided  a  basis  upon  which  such  a  theory  could  be 
proposed.  The  approach  adopted  by  both  G.  B.  Kistiakowsky  (5)  and  A.  R* 
Ubbelohde  ($)  is  that  of  a  preliminary  deflagration  evolving  discontinu- 
ously  into  a  detonation. 

Consider  a  mass  of  e:q)losive  particles  in  which  a  small 
peripherial  deflagration  has  been  started.  As  the  deflagration  spreads 
inward,  pressure  gradients  develop  due  to  the  restricted  flow  of  the 
gaseous  products  through  the  interstices  between  the  particles.  As  the 
interior  pressure  rises,  the  speed  of  the  deflagration  increases,  result- 
ing in  a  faster  decomposition  and  a  further  increase  in  pressure.  This 
constantly  increasing  deflagration  may  evolve  into  a  stationary  combus- 
tion in  which  the  flow  of  heat  and  materisJL  away  from  the  reaction  zaie 
compensates  for  that  produced  by  the  reactions,  or  the  deflagration  rate 
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may  continue  to  increase  until  all  the  material  is  consumed  provided  the 
mass  of  explosive  involved  is  less  than  the  "critical  size".  The  concept 
of  "critical  size"  for  a  particular  e;xplosive  is  a  functicai  of  many 
variables  includirig  shape  (volume-area  ratio),  particle  size,  crystal 
form  for  polymorphic  substances,  and  density.  It  is  probably  best 
defined,  although  somewhat  nebulously,  as  a  limiting  criterion,  that  is, 
the  critical  size  is  that  size  at  which  the  rate  of  decomposition  under- 
goes a  discontinuity.  The  discontinuity  may  result  in  a  partial  or 
complete  extinguishing  of  deflagration,  with  perhaps  a  re-initiation  at 
an  adjacent  spot  in  the  sxibstance  or  it  may  be  evidenced  by  the  phenomenon 
of  detonation.  To  see  how  the  non-stationary  deflagration  might  evolve 
into  a  detonation  we  must  consider  the  flow  of  gases  from  a  center  of 
deflagration  in  more  detail.  As  each  crystal  decomposes  it  acts  as  a 
center  from  which  the  gas  molecules  stream  outwards,  thus  forming  a 
pressure  wavelet.  If  a  number  of  these  wavelets  can  integrate  to  give 
a  streaming  velocity  corresponding  to  that  accompanying  a  stable 
detonation  wave  for  the  system,  we  have  conditicttis  identical  to  those 
observed  in  a  detonation.  The  type  of  flow  just  postulated  has  a  very 
close  analogy  in  the  field  of  gaseous  combustion.  A  column  of  gas  acted 
upon  by  a  piston  subject  to  acceleration  is  frequently  cited  as  eui  example 
of  the  mechanism  of  shock  wave  formation  in  gases.  Further  discussion  of 
this  mechanism  is  included  in  the  following  section. 

Assuming  a  shock  wave  has  developed  it  is  of  interest  to  con- 
sider how  its  passage  through  the  e::qplosive  can  initiate  the  decomposi- 
tion. The  details  of  the  mechanism  by  which  a  shock  with  its  discontin- 
uous rise  in  pressure  and  non-isentropic  temperature  effects  activate 
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the  explosive  molecules  is  still  the  subject  of  considerable  specxila- 
tion. 

It  has  been  proposed  by  Ubbelohde  (5)  that  the  most  probable 
mechanism  for  the  propagation  of  detonation  is  the  collision  of  high 
speed  atoms,  radicals  or  molecules  from  behind  the  vraive  front  with  the 
quiescent  layer  of  molecules  just  ahead  of  the  front.  Particles  behind 
the  detonation  front  are  assumed  to  have  an  unusually  high  activation 
efficiency  due  to  high  random  velocities  at  the  detonation  temperature 
and  a  high  directed  velocity  resulting  from  the  hydrodynamic  conditions 
of  propagation  of  a  shock  wave.  Impact  of  a  high  speed  paurticle  is 
visualized  as  disrupting  an  unreacted  molecule  before  the  kinetic  energy 
has  had  time  to  be  distributed  among  all  the  vibrational  degrees  of 
freedom.  If  dissociation  occurs  after  the  redistribution  of  energy, 
then  such  a  dissociation  must  be  attributed  to  high  temperature  in  the 
ordinary  sense.  Since  redistribution  requires  a  small  or  zero  activation 
energy,  whereas  dissociation  requires  an  activation  energy  about  equal 
to  the  strength  of  the  bond  broken,  it  is  probable  that  redistribution 
of  the  kinetic  energy  will  occur  before  dissociation  in  most  cases. 

It  has  been  suggested  by  Wendlandt  (6),  tyring  (7)  and  others 
that  the  mechanism  of  chemical  activation  depends  primarily  on  the  high 
temperature  produced  by  the  shock  wave.  Although  considerable 
uncertainty  exists  in  regard  to  reactions  at  the  extreme  temperatures 
and  pressures  of  detonation,  it  has  been  shown  (8)  that  the  temperature 
coefficient  of  the  energy-releasing  process  is  much  smaller  than  the 
usual  Arrhenius  value.  This,  plus  the  fact  that  it  is  xinlikely  that 
sufficient  time  exists  in  the  activation  process  for  the  vibrational 
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and  rotational  quanta  to  attain  the  full  Haxwell-Boltzman  equilibrium, 
indicates  that  the  temperature  effect  alone  probably  does  not  completely 
account  for  the  activation. 

Still  another  theory  has  been  proposed  suggesting  that  the 
pressure  in  the  wave  front  causes  a  bending  of  the  molecules  placing 
reactive  groups  in  close  proximity  to  each  other,  whereupon  the 
decompositicxi  occurs.  It  is  felt  that  this  theory  is  the  least  satis- 
factory of  those  presented  since  Bridgeman  (9)  has  shown  that  high 
hydrostatic  pressiire  alone  has  no  effect  on  explosives.  This  theory  is 
also  untenable  from  the  standpoint  of  chemical  mechanism  since  practically 
every  conceivable  first  step  in  the  decomposition  of  organic  nitro 
e^losives  is  endothermic  and  the  pressure  activation  concept  requires 
an  exothermic  initial  step. 

Thus  at  the  present  time  it  appears  most  reasonable  to 
attribute  the  propagation  of  detonation  to  the  energy  possessed  by  the 
shock  wave  as  a  whole  rather  than  to  attempt  a  more  detailed  and  perhaps 
erroneous  description  of  the  phenomenon. 
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SHOCK  WAVES 

The  theoretical  basis  for  the  development  of  a  theory  of 
detonation  stems  from  theory  of  shock  waves  proposed  by  Hugonlot  (10) 
In  the  nineteenth  century.  In  mathematical  Isuiguage  he  formulated  the 
difference  between  shock  waves,  or  compresslonal  waves  of  finite 
amplitude,  and  weak  acoustic  waves.  The  former  were  found  to  be 
propagated  with  velocities  which  are  greater  than  the  acoustic  velocity 
of  the  medium,  the  more  so  the  higher  the  peak  pressure  in  the  crest  of 
the  wave.  The  medium  itself  is  moving  in  the  shocks,  the  wave  velocity 
being  slightly  less  than  this  mass  velocity  plus  the  acoustic  velocity. 
Regardless  of  the  original  shape  of  the  (finite)  compresslonal  wave  in 
the  medium,  a  shock  front  is  formed  eventually  as  the  wave  propagates 
through  the  medium.  This  means  that  pressure,  density,  and  temperature 
suffer  discontinuous  increases  in  the  front  of  the  wave.  In  the  course 
of  its  propagation  through  a  medium  a  pure  shock  wave  decays,  this 
process  being  superimposed  upon  the  normal  weakening  of  waves  due  to  the 
geometric  factor  of  divergence,  unless  mechanical  energy  is  constantly 
fed  into  the  shock  wave.  In  other  words,  the  mechanical  energy  of  shock 
waves  is  being  degraded  and  is  then  found  as  the  heat  energy  of  the 
medium  after  the  passage  of  the  wave. 

Since  a  detonation  wave  is  basically  a  self-sustaining  shock 
wave,  a  brief  qualitative  discussion  of  shock  wave  formation  appears 
desirable.  It  is  well  established  that  in  a  compresslonal  wave  of 
finite  amplitude  the  velocity  of  a  portion  of  a  compresslonal  pulse  is 
directly  proportional  to  its  amplitude.  Those  portions  of  the  pulse 
with  the  greatest  amplitude  will  travel  faster  than  portions  of  lower 
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amplitude,  thereiore  changing  the  shape  of  the  pulse  as  indicated  in 
Figure  1. 


->  b 


Figure  1 
Formation  of  a  Shock  Front  Pb>Pa 

Becker's  treatment  (11)  of  the  origin  of  a  shock  vfave  in  a 
gas  is  also  of  interest.  Let  a  piston  in  a  long  tube  be  given  a  sudden 
small  velocity.  The  gas  in  front  of  the  piston  vdll  be  compressed  and 
a  pressure  wave  vdll  move  ahead  of  the  piston  with  the  velocity  of  sound 
in  the  medium.  The  gas  in  front  of  the  piston  and  behind  the  pressure 
wave  will  be  moving  with  the  velocity  ol  the  piston.  If  the  velocity 
of  the  piston  is  increased  by  another  small  increment,  a  new  pressure 
wave  arises  that  moves  somewhat  faster  than  the  first,  since  it  moves  in 
the  gas  already  heated  by  the  foregoing  compression  and  moving  vdth  the 
velocity  of  the  piston.  Jiach  succeeding  wave  will  move  more  rapidly 
than  the  preceding  one  and  will  ultimately  overtake  the  initial  one  to 
form  a  composite  wave  with  an  infinitely  steep  front;  i.e.,  a  discon- 
tinuity in  pressure  and  temperature  will  be  formed. 

The  qualitative  argiiment  can  be  made  more  explicit  by  the 
Riemann  formulation  of  the  equations  of  hydrodynamics.  Such  a  treatment 
is  presented  by  Lamb  (12)  and  the  results  confirm  the  conclusions  reached 
above • 

By  reversing  the  argument  given  for  compression  waves,  it  is 
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apparent  tnat  a  rarefaction  wave  of  finite  amplitude,  in  which  the  later 
portions  of  the  wave  are  regions  of  lower  pressure,  will  spread  out  as 
it  advances  and  a  discontinuity  will  not  arise » 

Returning  to  the  shock  front,  which  for  the  purposes  of  this 
discussion  can  be  treated  as  a  mathematical  plane,  it  becomes  desirable 
to  establish  equations  relating  the  conditions  on  both  sides  of  the 
discontinuity.  The  equations  for  the  conservation  of  mass,  momentum  and 
total  energy  have  been  employed  by  Rankine  (13)  and  Hugoniot  (10)  in  the 
formulation  of  three  conditions  relating  the  pressure,  specific  volume 
and  particle  velocity  beidnd  a  discontinuity  moving  with  a  velocity  D  to 
those  properties  in  the  undisturbed  medium  aJiead  of  the  discontinuity. 

Consider  a  shock  front  moving  into  a  region  in  which  the 
particle  velocity  is  u^,  the  density  Pq  and  the  pressure  P  from  a 
region  in  which  the  same  properties  are  denoted  by  Uj_,(>2^  and  P^. 


Pl»fl>  «] 


Po>  Po'  ^o 


shock  front 
Figure  2 

Conditions  at  a  Shock  Front  Moving  with  Velocity  D. 

In  a  time  dt  a  mass  of  fluid  ^Q(D-UQ)dt  will  enter  unit  cross-section 
area  of  the  front  and   a  mass  P]_(D-u-]_)dt  will  leave  the  front  (see 
Figure  2).  As  dt  is  made  very  small  so  that  the  layers  on  either  side 
are  infinitesimal,  the  mass  flow  into  the  front  must  equal  the  mass  flow 
away  from  it.  Thus  an  expression  for  the  conservation  of  mass  may  be 
written: 
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The  change  in  momentum  associated  with  the  mass  flow  across 
the  frait  must  equal  the  net  force  per  unit  area.  Thus  for  the  con- 
servation of  ffl«mentum: 

fo(D-Uo)(ui-u^)  =  Pi-P^  (2) 

Fina4.1y  the  conservation  of  energy  requires  that  the  net  work 
doneby  the  pressures  P^   and  ?q   equals  the  increase  in  kinetic  and 
potential  energy.  The  work  done  by  each  pressure  is  Pudt  for  unit  area, 
the  kinetic  energy  is  ^  u^  and  the  internal  energy  is  E  per  unit  mass. 
Therefore, 

Pl^l-Po^o  =  (©(l^-'^o)  lil-Mi(«f-^)  (3) 

The  above  equations  may  be  rewritten  in  the  form  below  by 
substitution  of  V  =  l/P  and  some  algebraic  manipulation. 

D  '  ^^0*^0  y  CPl-PoV(Vo-Vi)  (4) 

^1  =  M  yj   (Pi-PoKVq-Vi)  (5) 

El-Eo  =  ^  (Po-^Pl)(VVi)  (6) 

The  three  equations  above  are  often  referred  to  as  the  Rankine- 
hugoniot  conditions. 

Before  proceeding  to  the  consideration  of  detonation  processes, 
it  is  desirable  to  consider  Equation  6  in  more  detain  arid  to  investigate 
the  thermodynamic  properties  of  the  shock  wave.  Knowledge  of  the 
equation  of  state  and  the  heat  capacity  enables  one  to  calculate 
E^-Eq  as  a  function  of  P^  and  V^  so  that  Equation  6  becomes  the  Rankine- 
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Hugonlot  equation  relating  the  pressux*e  and  volume  in  a  compressional 
wave  of  finite  amplitude  for  given  initial  conditions.  The  relation 
describes  the  change  in  internal  energy  before  and  after  compression  by 
the  shock  wave  as  a  fiinction  of  the  pressures  and  specific  volumes. 
This  change  is  adiabatic,  since  no  heat  is  exchanged  with  the  surroundings 
but  it  is  not  isentropic.  The  non-isentropic  nature  of  a  shock  (or 
detonation)  wave  was  not  appreciated  by  many  of  the  early  investigators 
and  led  to  a  conflict  with  the  energy  principle,  A  method  of  determin- 
ing the  increase  in  entropy  occurring  across  a  shock  front  is  discussed 
along  with  other  detonation  zone  parameters  in  a  later  section. 
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DKTQNATIQN  WAVES 

The  existence  o£   detonation  waves  was  originally  reported  by 
Berthelot  and  Vieille  (14)  in  connection  with  studies  of  gaseous 
ejqplosions.  Such  a  wave  differs  from  a  shock  wave  in  two  respects.  In 
the  first  place,  the  wave  is  sustained  and  propagated  through  the 
mechanism  of  a  chemical  reaction  supplying  energy.  Therefore,  the 
chemical  energy  released  must  be  included  in  the  calculations  of  the 
Rankine-Hugoniot  conditions.  The  second  difference  is  that  the  velocity 
of  the  wave  is  determined  by  the  thermodynamic  state  of  the  material 
immediately  behind  the  wave  front  rather  than  by  the  motion  of  an 
external  boundary  surface  as  is  the  case  in  pure  shock  waves.  The  three 
conditions  relating  properties  of  the  material  on  either  side  of  the 
detonation  front  in  terms  of  the  conservation  of  mass,  momentum,  and 
energy  are  no  longer  sufficient  to  determine  the  unknown  quantities 
behind  the  front  in  addition  to  the  velocity  of  the  front.  The  problem 
differs  from  that  considered  in  the  previous  section  in  that  there  is 
no  controllable  external  condition,  such  as  the  velocity  of  a  piston, 
available  to  determine  one  of  the  variables  behind  the  front. 
The  Chapman~Jouguet  Condition 

An  additional  relation  is  required  to  make  possible  a  unique 
solution  of  the  hydrodynamic  equations.  Chapman  (15)  postulated  a 
relation  between  the  velocity  of  the  detonation  wave  and  the  sonic 
velocity  in  the  gas  behind  the  wave.  He  also  assumed  that  the  minimum 
detonation  velocity  consistent  with  the  Rankine-Hugoniot  conditions  pre- 
vailed in  a  steady  detonation.  Although  his  final  results  were  in  error 
due  to  assuming  that  isentropic  conditions  prevailed  at  a  discontinuity. 
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the  basic  concept  was  identical  to  the  mathematical  expression  published 
by   Jouguet  (16)  two  years  later.  The  condition  now  bears  the  name 
Chapman- J ouguet  and  postulates 


D  =  Ci+U;l 


(7) 


where  \x^   is  the  material  velocity  and  Cj^  the  sordc  velocity  in  the 
reaction  products  behind  the  wave  front. 

Although  a  valid  theoretical  proof  of  the  necessity  of  this 
condition  is  not  yet  available,  there  is  little  doubt  that  it  correctly 
describes  the  observed  phenomenon.  Several  plausible  arguments  have 
been  advanced  which  aid  in  the  acceptance  of  this  condition o  The  general 
form  of  the  pressure-volume  curve  representing  possible  states  of  the 
products  after  detonation  is  shown  in  Figure  3.  This  curve  is  specified 
by  the  Rankine-Hugoniot  conditions,  the  thermochemical  data  concerning 
the  e^losive  and  the  equation  of  state. 


R-H  curve  (PiV^) 


A  (Po,  V^) 


V 

Figure  3 

Generalized  Rankine-Hugoniot  Curve 


The  point  A  represents  the  initial  conditions  existing  in  the 
explosive  before  the  passage  of  the  detonation  wave.  If  the  initial 
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material  velocity  is  zero,  as  is  invariably  the  case  vdth  solid 
explosives,  Lqpiation  4  gives  the  detonation  velocity  in  terms  of  the 
initial  and  final  states. 

D  =  Vo  ;  (Pi-Po)/VVi)  (8) 

From  Figure  3  it  is  evident  that  this  relation  can  be  written  as 

D  =  Vo  J  tan  9  (9) 

where  9  is  the  angle  between  the  system  of  lines  passing  through  point  A 
and  the  (negative)  V  axis.  For  any  value  of  0  greater  than  0^,  there 
will  be  two  distinct  final  states  (C  and  G)  corresponding  to  the  same 
value  of  D.  For  0=0^,  the  line  AB  is  tangent  to  the  curve  at  point  B 
and  is  the  fineil  state  corresponding  to  the  minimum  value  of  D. 

If  the  final  state  of  the  products  corresponds  to  a  point  C 
on  the  R-H  curve  higtier  than  the  point  B,  it  can  be  shown  (see  Appendix  A) 
that  the  velocity  of  the  detonation  wave  is  less  than  u-t-c  and  a  rare- 
faction wave  travelling  at  sonic  velocity  with  respect  to  the  products 
would  move  faster  than  the  detonation  and  hence  engulf  aind  destroy  it. 
(a  rarefaction  wave  will  arise  from  any  process  that  abstracts  energy 
from  the  reaction  products  such  as  heat  loss  or  turbulence.) 

At  the  point  B,  the  velocity  of  the  detonation  wave  is  equal 
to  the  velocity  of  sound  in  the  reaction  products  plus  the  mass  velocity 
of  those  pzx>ducts,  thus  a  rarefaction  wave  proceeding  with  its  maximum 
velocity  will  be  unable  to  overtake  the  detonation  wave.  The  mathe- 
matical proof  that  D  =  c<«i  at  point  B  follows  that  given  by  Kistiakowsky 
and  Wilson  (17).  From  Equations  4  and  5  the  detonation  velocity  is  given 
by 
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The  velocity  of  sound  in  a  gas  is  given  by  the  illation  (18) 

Tvls 


c  -y  (ap/ap)s  «  v(-- 


(10) 


(11) 


If  there  is  a  point  on  the  Rankine-Hugoniot  (R-H)  curve 
having  the  same  slope  as  the  corresponding  point  on  the  isentropic  curve 
(PV  curve  with  constant  entropy),  then  at  that  point 

( avjr  '  1^)3 

where  the  subscript  r  refers  to  the  Rankine-Hugoniot  curve.  Differ- 
entiation of  the  Rankine-Hugoniot  relation ,  £quation  6,  holding  P^  and 
Vq  constant  yields 

dEi  =  -  i(Po-*Pl)^Vifi(V^-V^)dP^         (12) 

Now,  from  the  first  law  of  thermodynamics,  dE  «  TdS-PdV,  thus  Equation 
12  can  be  written 


Tl 


feL= 


i(Vo-Vi) 


/£i  .  ^o-Pi 


(13) 


It  follows  that  the  slope  of  the  R-H  curve  through  a  point  at  which 
dS^  s  0  can  be  expressed  as  shown  in  Equation  14  below 


/8Pl\    ^  ld?i 

l^vi/r     lavi 


V^i 


tan  9^         ^^) 


This  is  also  the  slope  of  the  line  AB  which  is  tangent  to  the  R-H  cui^e 
at  point  B.     Now  the  detonation  velocity  at  any  point  on  the  R-H  curve 
is  given  by  Equation  4  and  is 


(A) 
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Therefore,  at  point  B 

D  =  V^  f-J^]^  =  Zo      (by  Eq.  11)   (15) 

Combining  this  vdth  Equation  10  yields 

D  =  u^+Cj^ 

The  argument  which  is  used  to  exclude  points  on  the  R-H  curve 
below  B  is  based  on  the  entropy  of  the  products.  Consider  the  two  points 
C  euid  G  lying  on  the  same  straight  line  from  A  and  corresponding  to  the 
same  detonation  velocity  D.  The  entropy  of  products  at  C  is  greater  than 
at  Q  since  the  transition  from  G  to  C  is  equivalent  to  passage  from  an 
initial  state  G  across  a  shock  wave  of  velocity  D  to  the  final  state  C* 
Since  there  is  a  definite  increase  in  entropy  accompanying  the  passage 
of  a  shock  wave  the  entropy  at  C  will  be  greater  than  at  G,  therefore 
the  upper  point  will  be  a  thermodynamically  more  probable  state  than  G. 
This  is  true  for  all  angles  of  Q>  &i» 

Thus  points  above  B  aire  more  probable  than  points  below  but 
points  above  B  are  not  stable  due  to  the  effect  of  the  rarefaction  wave. 

The  question  of  the  validity  of  the  Chapman- Jouguet  condition 
has  also  been  discussed  by  von  Neumann  (19)*  By  applying  the  conserva- 
tion Equations  4>  5  &nd  6  not  only  to  the  C(xapletely  reacted  material 
but  also  to  each  stage  of  the  reaction,  he  has  been  able  to  substantiate 
the  validity  of  the  condition  in  certain  cases.  In  addition  to  the 
single  R-H  curve  in  Figure  3  representing  the  final  state  of  the 
products,  there  will  be  a  series  of  R-H  curves  corresponding  to  the 
intermediate  stages.  Since  the  detonation  velocity  is  the  velocity  with 
which  the  region  of  chemical  reaction  progress,  it  is  the  same  for  each 
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intermediate  stage.  Therefore,  the  line  from  the  initial  state  P^,  V^ 
to  each  P^^,  V^^^   representing  the  condition  of  a  given  layer  must  have 
the  same  slope  and  must  therefore  be  the  same  line  for  all  degrees  of 
completion  of  the  reaction.  Since  in  order  to  get  from  incomplete  to 
complete  reaction,  the  material  must  have  traversed  all  intermediate 
stages,  this  line  must  touch  all  of  the  R-H  curves  corresponding  to 
different  degrees  of  completion.  The  family  of  R-H  curves  corresponding 
to  the  degree  of  completion  of  the  chemical  reaction  from  n  =  0  to 

n  =  1  has  the  general  form  shown  in  Figure  4. 

Cl 

Reaction  Complete 


Figure  4 
Family  of  Rankine-Hugoniot  Curves  from  n  =  0  t-o  n 


=  1. 


The  motion  of  a  point  on  the  diagram  representing  the  behavior  of  an 
element  of  explosive  passing  through  the  shock  front  and  into  a  finite 
reaction  zone  must  always  remain  on  a  straight  line  to  satisfy  the 
conservation  equations.  In  the  shock  wave,  the  representative  point 
must  begin  on  the  curve  n  =  0,  at  point  A,  must  intersect  the  curves 
for  increasing  n  in  order  and  finally  end  somewhere  on  the  curve  n  =  1, 
In  order  that  there  exist  a  discontiniiity  in  pressure  and  temperature 
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at  the  shock  front,  there  must  be  another  intersection  of  the  straight 
line  vdth  the  curve  n  =  0,  corresponding  to  point  C.  Thus  in  the  shock 
front  the  representative  point  runs  from  A  to  the  limiting  value  at  the 
point  C  where  it  again  meets  the  curve  n  ■  0.  So  far  no  chemical 
reactions  have  occurred  but  the  eaqjlosive  material  has  been  "shocked" 
into  a  condition  where  reactions  can  proceed,  ^iow  the  representative 
point  must  proceed  continuously  from  the  curve  n  =  0  to  n  =  Ij  this 
requirement  immediately  restricts  the  possible  intersections  to  points 
on  n  •  1  at  or  above  B.  If  the  straight  line  intersects  the  curve 
n  >B  1  at  a  point  above  B,  the  tangent  to  the  curve  at  that  point  would 
be  steeper  than  at  B,  so  the  velocity  of  sound  in  the  products  would 
exceed  the  detonation  velocity,  and  a  rarefaction  wave  would  engulf  and 
destroy  the  detonation  wave.  Thus  the  intersection  with  the  curve 
n  =  1  must  be  at  B  giAd.ng  a  point  of  tangency  and  corresponding  to  Chapman- 
Jouguet  condition,  D  =  u^^-^-c^. 

Von  Neumann  £ilso  discusses  the  case  of  intersecting  R-H  curves, 
but  since  this  possibility  does  not  correspond  to  any  known  physical  sit- 
uation it  will  not  be  pursued  further. 

By  a  straightforward  thermodynamic  procedure,  Scorah  (20)  has 
shown  that  when  the  Chapnan-Jouguet  condition  is  assumed  the  ratio  of 
entropy  increase  to  available  energy  passing  through  the  shock  front  is 
a  maximum.  Although  this  result  is  in  accord  with  the  fundamental 
principle  of  maximum  degradation  of  energy  observed  in  natural  processes, 
it  is  regarded  only  as  a  further  substantiation  of  the  condition  rather 
than  a  rigorous  proof. 
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TH£  DETONATION  ZONE 

Before  proceeding  vdth  the  deyelopanent  of  the  fundamental 

theory,  it  is  advantageous  to  have  an  idea  of  the  order  of  magnitude  of 

the  dimensions  of  the  zones  in  which  the  process  of  detonation  occurs. 

The  essential  elements  of  a  steady  one-dimensional  detonation  process 

can  be  illustrated  as  in  Figure  5*  The  detonation  zone  is  considered 

to  be  stationary  in  space  and  the  e^losive  feeds  into  the  zone  from 

riglit  to  left  at  the  detonation  velocity  D. 

^ Det.  Zone   » 


Steady  Plane  Detonation  in  a  Solid  Explosive 


Conditions  throughout  the  zone  are  considered  to  be  steady.  The  material 
first  passes  into  a  plane  shock  zone  where  abrupt  Increases  of  pressure , 
temperature,  mass  velocity  and  entropy  occur.  This  transition  takes 
place  within  a  narrow  region  whose  length  is  determined  by  the  viscosity 
and  the  thermal  conductivity.  If  the  dimensions  of  this  shock  zone  are 
very  small  in  comparison  with  other  lengths  under  consideration,  the 
shock  may  be  treated  as  a  mathematical  plane  or  discontinuity  without 
error.  For  one-dimensional  motion,  the  Navier-Stokes  equations  (21)  of 
motion  and  energy  for  a  steady  state  in  which  the  shock  wave  is  at  rest 
are: 


V   V, 


^dx 
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and 

ZiE  -  AQ  -i(D-u)2  -Po(Vo-V)  =  k  I2  ^ 

"  ^       D  dx 

where  /ny  is  the  viscosity,  k  is  the  thermal  conductivity,  u  is  the 
material  velocity,  T  is  the  temperature  and  x  is  the  coordinate  in  the 
direction  perpendicular  to  the  shock  front.  These  differential  equations 
have  not  as  yet  been  solved  in  closed  form  for  general  values  of  the 
parameters.  For  the  case  of  a  pure  shock  wave  (  ^Q  ■  0)  in  liquid  di- 
ethyl ether,  Becker  (22)  has  calculated  the  length  of  the  zone  to  be 
6.5A^  at  a  shock  pressure  of  10, OCX)  atms.  and  1,UA?   at  a  pressure  of 
100,000  atms.  For  a  solid  e:}qplosive  with  an  assumed  viscosity  of  0.2 
centipose  and  a  thermal  conductivity  of  10^  ergs  cm"  sec  ,  i^yring  (23) 
calculated  the  length  of  the  zone  to  be  of  molecular  dimensions.  Calcu- 
lations including  the  effects  of  viscosity  and  thermal  conductivity  in 
the  conservation  equations  have  shown  these  effects  to  be  of  importance 
only  at  the  very  front  of  the  detonation  wave,  where  the  amount  of  chem- 
ical reaction  is  negligible. 

The  distance  between  the  points  of  initiation  of  the  chemical 
reaction  and  the  point  at  which  the  Chapman- J ouguet  condition  is 
applied  is  called  the  reaction  zone.  There  are  a  number  of  methods  for 
estimating  the  length  of  this  zone.  Since  each  method  rests  on  the 
interpretation  of  a  different  phenomenon,  agreement  among  results  found 
by  different  methods  greatly  enhances  the  accuracy  of  the  determination. 
The  results  of  (1)  direct  observation  of  the  luminous  zone,  (2)  extrapo- 
lation of  low-temperature  rate  measurements,  (3)  surface  reaction  rate 
calculations  for  granular  explosives,  and  (4)  studies  of  the  building-up 
of  a  low  velocity  wave  to  steady  state  conditions,  all  indicate  that  the 
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reaction  zone  for  TNT  under  steady  state  conditions  is  about  ^  centimeter 
in  length.  Although  the  exact  figure  vdll  vary  from  explosive  to  explo- 
sive, it  is  reasonably  well  established  that  the  reaction  zone  length 
for  solid  ejq>losive3  is  between  1  mm.  and  1.5  cm. 

The  time  required  for  a  particle  to  pass  through  the  reaction 
zone  can  be  calculated  knowing  the  zone  length,  detonation  velocity,  and 
the  material  velocity.  Reaction  times  of  the  order  of  magnitude  of 
1  microsecond  have  been  found  for  practically  all  organic  explosives. 
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THE  HYDRQDYNAMIC  THEORY 

Application  of  the  Rankine-Hugoniot  relations  and  the 
Chapman- Jouguet  condition  to  an  actual  chemically  supported  detonation 
permits  evaluation  of  the  energy  relation.  Equation  6,   in  terms  of  heat 
of  reaction  and  specific  heat  data  for  the  reaction  products.  Since  the 
shock  zone  is  of  molecular  dimensions  any  contribution  of  heat  by  thermal 
conduction  is  nei 


jglected.  The  energy  relation  can  now  be  written 

To    *»     ^1 
or  introducing  the  mean  specific  heat  at  constant  volume  and  assuming  an 


ideal  gas 

i^ere  Q  is  the  heat  absorbed  in  the  reaction  of  M  grams  at  the  initial 
temperature  Tq*  As  indicated  above  the  computation  of  conditions  about 
the  detonation  front  requires  the  knowledge  of  the  following: 

1.  Equation  of  state  of  the  products 

2.  Composition  of  the  reaction  products 

3.  Heat  capacity  of  the  products 
4*  Heat  of  reaction 

In  the  case  of  solid  explosives,  the  selection  of  an  equation  of  state 
for  products  at  temperatures  of  the  order  of  3000  -  5000°  A  and  at 
pressures  of  50,000  to  200,000  atmospheres  presents  a  formidable  problem. 
Our  present  day  knowledge  of  equilibrium  constants  under  the  above 
conditions  is  also  quite  meager.  In  view  of  this  situation,  it  appears 
desirable  to  test  the  validity  of  theory  under  the  less  stringent  con- 
ditions provided  by  gaseous  explosions. 
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Gaseous  Kxploslves 

The  approach  adopted  in  this  discussion  is  similar  to  that 
applied  to  mixtures  of  hydrogen  sind  o:<ygen  by  Lewis  and  Friauf  (24). 
Since  the  equations  used  for  idesd  gases  will  later  prove  of  value  in 
the  preliminary  treatment  of  solid  explosives,  they  will  be  fully 
developed  here.  The  Rankine-Hugoniot  conditions  for  conservation  of 
mass  and  momentum,  Equations  1  and  2,  may  be  written  as  below,  if 
\Xq   is  assumed  to  be  zero. 

0  =  kHH  (18) 

Vo   Vi 


Dui  . 


=  Pn-Pn  (19) 

o 


Vr      1   ° 


KLimination  of  u^  from  the  above  equations  yields 

v: 


VPo  =  Z2   ( V^l)  (20) 


o 

At  the  Chapman-Jouguet  point  on  the  R-H  curve  the  detonation  velocity 
expression  Equation  4  can  be  written 

D  =  Vo  V  (Pi-PoV(Vo-Vi)  =  V^  7  -(clPi/dVi)3       (21) 

For  an  ideal  gas,  the  isentropic  e^ansion  law  is 

Pj^Vjl        =  constant 

Therefore,  differentiating  and  substituting  into  Equation  21 

D  =  Vo/Vi   ^nPiVi  =  Vq/Vi  JnniRTi/M        (22) 

in  which  n^^  is  the  number  of  moles  of  gas  per  M  grams  of  product  gases < 
Since  the  initial  pressure  P^  is  negligible  in  comparison  with  P^  in 
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Equation  20,  substitution  of  Equation  20  into  Equation  22  with  the 
elimination  of  D  yields 

X  =  FT-  ^'^^ 

^o  ^1 
This  expression  may  now  be  substituted  into  the  energy  relation. 
Equation  17,  again  neglecting  P  to  give 

If  the  chemical  reaction  under  consideration  went  quantitatively  and 
the  distribution  of  at<xiis  in  the  products  was  independent  of  temperature 
and  pressure.  Equation  24  could  be  solved  for  the  final  temperature  T!-^ 
if  the  variation  of  Yj^  and  Cy  with  T  were  known,  Actvially,  of  course, 
both  n^  and  Q  depend  upon  pressure  and  temperature.  The  calculation  is 
therefore  somewhat  more  tedious  but  the  principle  involved  is  unchanged. 

Once  the  final  temperature  is  found,  the  detonation  velocity 
can  be  computed  from  the  combination  of  Equations  22  and   23 

D  =  (Yi+l)  y  niRTi/riM  (25) 

Lewis  and  Friauf  have  calculated  detonation  velocities  for  mixtures  of 
hydrogen  and  oxygen  with  various  excesses  of  both  components  and  with 
the  addition  of  inert  gas.  They  allowed  for  the  equilibria 
^  H2+O2  =  H2O,  Hg  =  2H,  and  H20f^02  =  20H.  Although  the  dissociation  of 
O2  was  not  considered  and  no  allowance  was  made  for  the  excess  heat 
capacity  of  oxygen  molecules  due  to  electronic  excitation,  there  is  no 
reason  to  suggest  that  a  similar  calculation  undertaken  today  would 
differ  significantly  from  their  results.  The  following  table  shows 
their  calculated  values  and  the  experimental  results  obtained  by  Dixon 
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(25)  and  Payman  and  Walls  (26). 

Table  I 


Mixture 


Detonation 

Velocity 

D,  m/i 

sec 

Ccffliputed 

Observed 

Deviation  % 

2806 

2819 

-0.4 

2302 

2314 

-0.7 

1925 

1922 

0.2 

1732 

1700 

2.0 

2378 

2A07 

-1.2 

2033 

2055 

-1.1 

1850 

1822 

1.1 

3354 

3273 

2.5 

3627 

3527 

2.6 

3749 

3532 

6.1 

(2H2+02) 

(2H2K)2)  +  IO2 
(2H2+O2)  -k  3O2 
(2H2-K)2)  -f  5O2 
(2H2+O2)  4  IN2 
(2H2+O2)  +  3N2 
(2H2+O2)  +  5N2 
(2H2'K)2)  +  2H2 
(2H2+O2)  +  4H2 
(2H2+O2)  f  6H2 


The  agreement  between  the  cadculations  auid  the  observed  measurements 
serves  as  an  extremely  good  verification  of  the  hydrodynamic-thermo- 
dynamic  theory.  With  one  exception  the  differences  between  the 
comparative  values  are  all  within  the  experimental  errors  of  velocity 
measurements.  The  detonation  velocity  of  a  number  of  gaseous  ejqplosives 
of  the  hydrocarbon-oxygen  type  also  have  been  calculated  and  found  to 
be  in  good  agreement  with  experimentally  determined  values. 
Solid  Explosives 

The  extension  of  the  fundamental  theory  to  the  detonation  of 
solid  e:;qplosives  requires  the  selection  of  an  equation  of  state  applicable 
to  gases  under  extreme  conditions  of  temperature  and  pressure.  Use  of 
the  ideal  gas  law  is  immediately  ruled  out  since  the  experimentally 
observed  detonation  velocities  are  quite  sensitive  to  initial  density 
and  Equation  25  based  on  the  ideal  gas  law  shows  no  dependence  on  this 
factor. 

Another  complicating  factor  arising  with  solid  ejqplosives  is 
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the  question  of  the  composition  of  the  reaction  products.  To  consider 
all  of  the  various  equilibria  possible  when  a  typical  explosive  contain- 
ing atoms  of  carbon,  hydrogen,  nitrogen  and  oxygen  detonates  is  a 
tremendous  mathematical  task  and  cannot  be  justified  in  view  of  the 
uncertainty  of  the  thermodynamic  functions  of  the  possible  products. 
Justification  for  ignoring  the  more  complicated  equilibria  involving 
free  atoms,  polyatomic  molecules  and  free  carbon  in  approximate 
computations  will  be  presented  later. 
Variation  of  Parameters  -  Abel  Equation  of  State 

To  illustrate  the  form  assumed  by  the  hydrodynamic  equations 
and  to  ascertain  certain  qualitative  aspects  of  the  process  of  detona- 
tion a  particularly  simple  equation  of  state  will  be  chosen  as  an 
example.  An  equation  suitable  for  this  purpose,  but  of  less  value  in 
connection  vdth  practical  calculations,  is  the  Abel  equation, 
P(V-  «^)  =  RT,  with  the  covolume  o<   assumed  constant. 

The  details  of  deriving  the  relations  will  for  the  most  part 
be  omitted.  The  derivations  are  quite  straightforward  and  similar  to 
those  given  by  L^yring  and  his  associates  (7). 

Evaluation  of  the  sonic  velocity  appearing  in  the  Chapman- 
Jouguet  conditions.  Equation  7|  in  terms  of  the  equation  of  state  yields 

(D-U3^)2  .  Vi^y  _^1 (26) 

Vi-o< 

which  when  inserted  into  the  Rankin e-Hugoniot  relations  and  combinations 
of  these  relations  gives  the  following  five  equations.  (The  initial 
pressure  Pq  is  negligible  in  comparison  to  P;j^  and  is  omitted.) 
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(27) 


Pl(V  -o<)  ■  2R  ^"^yo  (28) 


Cy 


Q-fCvTo  2 


Gv    ^+1 


(29) 


^2  , 


D''  =  (Q+CyTo)  2(Y^^1)  ^  (30) 


(l-»-°^Vo)^ 


(31) 


Equations  27  to  31  furnish  the  detonation  velocity  and  all  other 
properties  of  the  reaction  products  provided  we  have  data  concerning 
the  heat  or  reaction,  the  mean  heat  capacity  of  the  products,  the  co- 
volume  and  the  gas  constant  per  gram.  It  is  apparent  that  neither  the 
detonation  velocity  nor  any  of  the  properties  of  the  products  depend 
upon  the  rate  of  chemical  reaction. 

An  estimate  of  the  probable  accuracy  of  the  calculated 
properties  can  be  made  from  the  above  equations.  The  final  temperature 
is  a  function  of  the  heat  of  reaction  and  the  mean  heat  capacity  of  the 
products  and  in  particular  is  independent  of  the  covolumej  thus  the 
temperature  should  be  determinable  to  a  hi^er  accuracy  than  properties 
depending  on  additional  data.  In  particular  the  final  pressure  is  least 
accurately  obtained  since,  in  addition  to  dependence  upon  the  thermal 
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data  and  the  covolumej  it  is  also  a  function  of  the  gas  constant  which 
depends  upon  the  mean  molecular  weight  of  the  products. 

The  properties  of  the  reacting  explosive  can  be  determined 
at  various  intermediate  stages  by  replacing  Q  in  Equations  2?  to  31  by 
nQ  where  n  denotes  the  degree  of  completion  of  the  reaction.  It  is 
interesting  to  note  that  such  a  procedure  is  another  application  of  the 
mechanism  used  by  von  Neumann  in  his  attempt  to  justify  the  selection  of 
the  Chapman-Jouguet  condition  based  upon  the  family  of  RanKine-Hugoniot 
curves.  The  variation  of  the  properties,  temperature,  pressure,  ratio 
of  c+u  to  D  and  entropy  of  a  typical  explosive  as  the  reaction  proceeds 
is  indicated  in  Figure  6.  BAch  property  is  plotted  in  terms  of  percent 
of  the  final  value. 
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Figure  6 
Variation  of  Properties  within  the  Detonation  Zone. 
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At  the  shock  frc»it  the  pressure  is  twice  the  value  attained  at  the  end 
of  the  reaction  zone,  since  the  substitution  of  nQ  into  Equation  28 
and  simplifications  yields,  to  close  approximation,  the  relation 

Pn  -  (1  +  /T^)Pi 

The  mass  velocity  yx^   is  always  the  same  fraction  of  its  final 
value  as  is  the  pressure  Fj^j  therefore,  the  same  curve  represents  the 
variation  of  both  parameters  through  the  reaction  zone. 

The  ratio  of  (c+u)  to  D  exceeds  2  at  the  shock  front  and  reaches 
unity  at  the  end  of  the  reaction  zoie,  thus  it  is  seen  that  within  the 
reaction  zone  a  rarefaction  wave  of  velocity  ou  will  continually  over- 
take and  decrease  the  velocity  of  a  discontinuity  progressing  at  a  lower 
velocity.  At  the  end  of  the  reaction  zone  D  =  c-^u,  (the  Chapman-Jouguet 
condition)  and  a  stable  detonation  wave  is  possible. 

The  temperature  and  entropy  show  sharp  increases  at  the  shock 
front  and  a  gradual  rise  throughout  the  reaction  zone. 
Equation  of  State  Based  on  Enerj^y  Considerations 

In  order  to  improve  the  agreement  between  calculated  and 
observed  values  of  the  detonation  velocity,  equations  of  state  more 
representative  of  the  extreme  conditions  confronted  have  been  suggested. 
Jones  (27)  has  developed  an  equation  suitable  for  use  at  high  pressures 
based  upon  Bridgman's  (28)  results  with  nitrogen  at  pressures  up  to 
15,000  atmospheres.  The  equation  takes  the  form 

where  EqCVn)  is  assumed  to  be  the  potential  energy  of  interaction  of 
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molecules  in  the  gaseous  product  and  N  is  the  number  of  moles.  The 
internal  energy  E^  is  expressed  in  the  form 

in  which  the  first  term  represents  the  repulsive  forces  and  the  second 
the  attractive  forces.  The  constant  C  is  assigned  a  value  to  make  the 
minimum  energy  equal  to  the  heat  of  vaporization.  The  function  f  is 
evaluated  from  compressibility  data  and  the  constants  A,  B  and  •<  from 
the  energy  of  the  normal  modes  of  vibration  of  molecules  in  the  fully 
excited  solid  state. 

At  lower  temperatures,  Jones  assumes  an  equation  of  the  virial 
form,  the  constants  being  adjusted  so  that  the  pressure  computed  at  an 
intermediate  value  of  V  by  the  both  relations  is  the  same.  Extending 
the  interaction  effects  determined  by  Bridgeman  to  higher  temperatures 
£uid  densities  the  total  energy  of  the  reaction  products  can  be  expressed 
ae 

E  =  Z[NiEi  +  (ZNi-N3)(Eo^-3/2RT)| 

where  Nj^  is  the  number  of  moles  of  each  molecular  species,  both  gas  and 
solid,  Ej^  is  the  energy  per  mole  exclusive  of  the  interaction  of 
molecules  in  the  gas  phase  and  N.  is  the  number  of  moles  of  solid 
products.  The  last  term,  therefore,  accounts  for  the  energy  of  molecular 
interaction  plus  the  vibration  energy  of  3/2RT  per  mole.  The  energies 
E^  can  be  evaluated  from  specific  heat  data  at  ordinary  pressures  since 
the  interaction  effects  at  high  pressures  are  represented  by  E^.  To 
determine  the  composition  of  the  reaction  products  as  a  function  of 
temperature,  the  activities  of  the  various  products  are  introduced. 
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The  activity  of  each  species  can  be  expressed  in  terms  of  the  composition 
N^  and  the  equilibrium  constants  for  the  assumed  reactions  as  a  function 
of  temperature  at  normal  pressures  to  which  ideal  gas  conditions  are 
applicable.  The  increase  in  internal  energy  of  the  products  must  be  the 
sum  of  the  chemical  energy  Q  released  by  the  reaction  and  the  work  done 
by  the  pressure  at  the  detonation  front.  Thus  a  knowledge  of  £,  Q  and 
the  equation  of  state  provides  sufficient  information  to  solve  the 
detonation  front  equations. 

Although  the  above  procedure  involves  several  major  assumptions 
to  make  it  feasible,  the  calculation  is  significant  because  the  results 
tabulated  in  Table  III  give  an  idea  of  the  success  obtainable  using  an 
equation  of  state  fitted  to  data  other  than  that  obtained  by  detonation 
measurements. 
The  Halford-Wilson-Kistiakowsky  Equation  of  State 

Information  on  gas  states  at  the  pressures  and  tempera txires 
encountered  in  detonations  can  be  obtained  by  the  use  of  the  hydro- 
dynamic  theory  if  the  variation  of  detonation  velocity  is  known  as  a 
function  of  the  initial  specific  volume  (loading  density) .  Since  the 
detonation  velocities  of  solid  ejqplosives  show  a  definite  increase  with 
loading  density,  such  measurements  provide  a  suitable  means  for  finding 
the  deviations  of  e:)qplosion  gases  from  ideal  gas  behavior.  The  detonation 
velocity-loading  density  relationship  is  of  little  value  in  suggesting 
the  correct  form  of  an  equation  of  state,  but,  once  a  form  is  established, 
the  relationship  is  sufficient  for  the  quantitative  evaluation  of  its 
empirical  parameters.  Unfortunately,  the  pressure-volume  relationship 
is  rather  insensitive  to  the  form  chosen  which  suggests  the  measurement 
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of  detonation  temperature  for  determining  the  appropriateness  of  the 
choice.  Since  detonation  temperatures  have  not  as  yet  been  determined 
to  a  degree  of  accuracy  sufficient  to  provide  a  check,  we  are  forced  to 
accept  the  form  which  fits  the  data  best  in  the  largest  number  of  cases. 

An  equation  of  state  found  to  be  very  successful  in  the  cal- 
culation of  detonation  velocities  is  one  developed  by  R.  S.  Halford  in 
cooperation  with  Wilson  and  Kistiakowsky  (17) •  The  equation  of  state 
has  the  following  form  when  all  the  products  of  reaction  are  gaseous 

PVM  =  nRT(l4xel^)  (32) 

with  X  «  Inik^/T^^VM 

V  is  the  specific  volume,  n  is  the  number  of  moles  of  gas  per  M  grams; 
P,  T  and  R  are  the  pressure,  temperature  and  the  gas  constant  per  gram 
mole,  respectively;  and  nj^  and  k^^   are  the  number  of  moles  and  an  empirical 
"covolume"  parameter  of  the  i-th  species,  respectively;  the  summation 
extends  over  all  gas  species  present.  Satisfactory  results  have  been 
obtained  with  the  constants  oi«  0.25  and  ^  =  0.30. 

In  cases  in  which  solid  proaucts  are  formed  in  the  detonation, 
the  equation  is  slightly  modified.  It  is  assumed  that  the  volumes  of 
gas  and  solid  are  adaitive  and  that  the  gas  obeys  Equation  32  independent 
of  the  presence  of  solid  components.  Since  thermal  expansion  and  baric 
compression  of  solids  have  opposing  effects  in  respect  to  the  volume,  and 
since  these  effects  are  hot  known  at  the  elevated  temperatures  and 
pressures  encountered,  it  is  assumed  that  the  specific  volume  of  a  solid 
component  is  independent  of  temperature  and  pressure.  Thus  the  specific 
volume  of  solid  components  may  be  assigned  values  determined  under  normal 
conditions  or  may  be  treated  as  adjustable  parameters.  With  these 
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assumptions,  the  equation  of  state  for  a  mixture  of  solid  and  gaseous 
components  becomes 

PM(V-Vs)  "  nRTF(x)  (33) 

where  V^  is  the  volume  of  solid  constituents  per  gram  of  mixture, 
obtained  by 

^s  =  na  ^a  "^  "^b  ^b  +  •••  ^^^*^ 

^  ,  rry^,.,.   are  the  weight  fractions  and  V^^,  V^,,..,  are  the  specific 
volumes  of  solid  constituents  a,  b,...  F(x)  is  a  function  defined  by 

F(x)  s  1  +  xe  ^^  (33b) 

with 

X  =  k/T*  M(V-V3)  (33c) 

k  =  rnik^  (33d) 

The  value  of  ec  and  ^  and  all  other  symbols  are  the  same  as  those  used 
with  Equation  32.  A  theoretical  evalxiation  of  the  covolume  factors 
(c.f.  Equation  33d)  is  impractical  in  view  of  the  extremely  high  density 
encountered  in  the  reaction  zone.  These  factors  are  therefore  determined 
by  comparison  of  the  results  of  a  calculation  of  detonation  velocity 
based  on  an  assumed  value  with  the  actual  observed  velocity.  As  an 
initied  approxLniation  only  the  most  abvuidant  detonation  products  are 
considered  and  the  total  covolume  constant  k  of  the  mixture  is  obtained 
for  several  explosives  at  a  number  of  loading  densities.  To  evaluate 
the  specific  covolume  factors  k^,  it  is  assumed  that  k  is  an  additive 
function  of  specific  values.  The  individual  values  are  then  adjusted 
by  an  iterative  procedure  for  optimum  conformance  with  the  mixture  values. 
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Since  the  original  work  was  completed,  other  methods  of  evaluating  these 
constants  have  been  suggested  and  a  wider  range  of  explosives  has  been 
used  as  a  basis,  yielding  slightly  modified  values.  It  has  been  found 
that  the  calculation  of  product  composition  is  much  more  sensitive  to 
the  individual  values  of  kj^  than  is  detonation  velocity. 

The  application  of  the  Halford-Wilson-Kistiakowsky  equation 
of  state  to  the  cedculation  of  detonation  velocities  and  related 
quantities  will  be  illustrated  below. 


CALCULATION  OF  THE  DETONATION  yJO^OCITY  AUD  REUTED  QUANTITIES 

The  calculation  of  the  detonation  velocity  and  related 
quantities  using  the  Halford-Wilson-Kistiaicowsky  equation  of  state  takes 
into  account  the  effects  of  gas  imperfections  and  results  in  an  entirely 
acceptable  set  of  parameters  in  the  vast  majority  of  explosives  to  which 
it  has  been  applied. 
Calculation  of  Detonation  Velocity 

From  pure  thermodynamics,  we  obtain  the  relation 

Differentiating  the  equation  of  state,  Equation  33,  holding  T 
constant  gives  the  result 

|<iP\   s  -  P(l»2xe  ^^4  (3x^e^^)  (-.rs 

From  Equation  33b,  the  coefficient  of  (V-Vg)  can  be  represented 

by  F(x)  and  that  of  P  by  y(x)  where 

y(x)  =  F(x)  f  xdF(x)  (36) 

dx 

Thus  Equation  35  can  be  written 

(dv)^  ■  (V^Vgi  F(x)  '"rldvjg  ^^'^^ 

Substitution  into  the  basic  equation  for  D,  Equation  21,  yields  the 
expression 

D  =  Vq  i  Pl^l  y(xi)/  (Vj^-Vg)  F(xi)         (38) 
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Introducing  the  equation  of  state  once  again  gives 

D=-!2^   yriniRTiy(xi)/M  (39) 

Equation  38  may  be  substituted  into  Equation  U  to  obtain  an  e:Q>re8sion 
in  terms  of  the  specific  volumes  (Pq  is  neglected) 

p  ,  Piy(xi)  yi(Vo-yi)  (40) 


which  on  rearrangement  becomes 

V^s   y(^)^i 


(U) 


Calculation  of  the  Specific  Heat  Ratio 

From  thermodynamics  J  one  obtains  the  follovd.ng  equations  for 
the  specific  heat  at  constaurit  volume  and  the  ratio  of  the  specific 
heats  at  constant  pressure  to  that  at  constant  volume 

C^  =  C*  +  m£'^s   (cl2p/dT^)y  dV  (42) 

where  Cy  is  the  ideal  heat  capacity  per  H  grams  and 


Y=  1 


-  'i  (f )!  /[%{  (^3) 


Introducing  the  equation  of  state  in  these  egressions  yields 

Cy  =  C*  +  nRo((2(x)  -  1)  (44) 

and  r  =  1  +  nR  [ z'(3E)2/c^(x)]  (45) 

where     z(x)  =  F(x)  -txx  ^F(xl  (/^) 

dx 
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Calculation  of  the  Temperature 

For  an  imperfect  gas,  Equation  l6  egresses  the  application  of 
the  first  law  of  thermodynamics 

M(El  -  Eo)  =  Q^  (Ti-T  )+M  (^^l"^8\cl£/dV)TT  *^^    ^^^^ 

(The  reaction  is  first  carried  out  to  the  final  temperature  at  a  large 
volume,  such  that  the  product  gases  are  ideal,  and  these  gases  are  then 
compressed  to  their  final  volume  vdth  a  resultant  energy  term  due  to  the 
gas  imperfections.) 

Introducing  the  equation  of  state  and  standard  thermodynamic 
relations,  one  obtains 

M  j^^l"^8(dE/dV)x^  dV  «  n;L^Tio<(F(x-L)-l)  (A8) 

Equation  /*!  can  be  rearranged  to  the  form 

V,.V^  =  ^(\)  CVi-Ve)  (49) 

which  when  combined  with  the  Hugoniot  Equation  6  yields 

Ei-Eq  =  ^PiF(x^)Vi/y(x;L)  y;L  =  ^niRTiF2(xi)/y(xi)T3_M  (50) 

Therefore,  from  Equations  47,  46  and  50,  one  obtains 

Ti  =  Cj  To-Q/C^+niR'^(F(xi).l)-n^RF2(xi)/2y(xi)y3^  (51) 

Solution  of  the  above  equations  is  based  on  an  assumed  knowledge  of  the 
final  composition,  thus  an  iterative  approach  appears  desirable.  If  a 
large  number  of  ejq^losives  are  to  be  investigated  a  two-stage  method  is 
recommended.  This  involves  the  calculation  of  ideal  values  of  D,  T,  and 
y  based  on  the  assumption  that  the  product  gases  are  ideal  and  then 
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evaluating  the  ratio  of  actual  to  ideal  values  using  tables  prepared 
for  this  purpose. 

The  entropy  increase  occurring  in  the  detonation  zone   can  be 
estimated  if  values  for  the  increase  in  temperature  and  volume  across 
the  shock  wave  are  assumed. 

The  entropy  increase  across  the  shock  zone  due  to  an  irreversible 
compression  to  a  pressxire  P  and  temperature  T  is  given  by 

The  entropy  change  due  to  the  irreversible  chemical  reacticai 
can  be  cailculated  by  proceeding  from  the  initial  to  the  final  state  by 
hypothetically  reversible  steps,  applying  the  relation 

^S  =  SVod.  -  S^i^los.  +  nCpln  Zi  -  n  R  In  ^1    (53) 

io        '^o 

in  iNtiich  the  subscript  1  refers  to  the  conditions  prevailing  at  the  end 
of  the  reaction  zone. 
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COMPOSITION  OF  THIi.  DETONATION  PRODUCTS 

In  order  to  carry  out  calculations  of  the  detonation  velocity 
and  other  pareimeters  of  interest,  it  is  necessary  either  to  calculate 
the  equilibrium  composition  of  the  reaction  products  or  to  make  an 
assumption  regarding  this  composition.  Before  discussing  the  detailed 
methods  used  to  determine  the  composition,  a  few  general  remarks  concern- 
ing equilibria  among  possible  products  appear  desirable* 

Dissociation  into  free  atoms  and  radicals  is  fortunately  not 
of  importance  at  the  high  pressures  encountered  in  detonation  of  solid 
ejqjlosions.  The  dissociation  constant  of  H2  =  2H  at  5000°K  has  been 
calculated  statistically  to  be  44*7  atm.  Since  the  concentration  of 
hydrogen  gas  rarely  exceeds  10^  by  volume  in  most  e^losives,  it  is  found 
that  at  a  pressure  10^  atmospheres  about  2%  of  the  hydrogen  Is 
dissociated  into  atoms.  This  corresponds  to  an  absorption  of  heat  of 
about  3  kilocalozdes  per  kilogram  of  explosive,  whose  total  heat  of 
explosion  is  of  the  order  of  1000  kcal.  per  kg.  Thus,  the  final 
temperature  is  lowered  to  the  extent  of  0.3^  by  the  inclusion  of  this 
dissociation  but  at  the  same  time,  the  number  of  moles  of  products  and 
the  pressure  is  increased.  The  two  effects  oppose  each  other  as  far  as 
detonation  velocity  is  concerned  and  it  is  estimated  that  the  net  effect 
results  in  a  0.10^  decrease  in  the  velocity.  The  dissociations  of  O2 
and  N2  occur  to  a  much  lesser  extent  and  can  also  be  ignored  in  all 
calculations.  Although  the  dissociations  of  water:  HgO  ■  ^  H2  -^  OH 
and  H^O  =  ^  O2  +  H2  take  place  to  only  a  very  small  extent  these 
equilibria  are  considered  in  precise  calculations. 

Various  polyatomic  molecules  have  been  detected  in  gas  samples    ^ 
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analyzed  after  closed  bomb  detonations,  Schmidt  (29)  has  reported  small 
but  variable  amounts  of  CH- ,  WHq,  C2H2*  HCN  and  C2N2  from  almost  all 
explosives.  It  is  suspected  that  the  majority  of  these  molecules  are 
formed  after  the  detonation  is  complete  and  while  the  products  are 
cooling.  Allowance  is  made  in  precise  computations  for  the  formation  of 
methane  and  ammonia,  but  the  remaining  three  are  neglected  in  the 
interest  of  mathematical  simplification.  If  it  is  assumed  that  the 
analytical  data  represents  the  true  conditions  existing  in  the  detonation 
wave^  the  formation  of  all  these  complex  molecules  has  a  small  but  not 
negligible  effect  on  the  results  of  calculation.  For  example,  in  TNT 
where  the  oxygen  deficiency  is  large  and  the  polyatomic  molecules  listed 
above  are  particularly  abundant  neglecting  their  formation  results  in  an 
overestimation  of  the  detonation  velocity  by  about  3.5^» 

The  formation  of  free  hydrogen  in  absence  of  solid  carbon  is 
considered  in  precise  computations  even  though  the  equilibrium  constant 
of  the  "water-gas"  reaction,  H2O  -«'  CO  =  C02+H2  ranges  from  0.2  at  3000° 
to  0.07  at  5000*^K.  Since  the  number  of  moles  of  products  does  not  change 
the  only  effect  on  the  calculations  is  that  due  to  the  evolution  of  heat. 
The  neglect  of  this  equilibria  in  the  absence  of  solid  carbon  can  result 
in  an  underestimation  of  the  rate  by  only  a  few  percent.  If,  on  the 
other  hand,  free  carbon  is  present,  the  reaction  C+H2O  =  CO+H2  proceeds 
almost  quantitatively.  The  equilibrium  constant  of  this  reaction  is  of 
the  order  of  10^  atm.  at  3000°  and  10^  atm.  at  5000°K.  Since  the 
fugacity  of  CO  in  the  reaction  zone  is  usually  less  than  10^  atm.,  the 
ratio  ^n^/^HjO    is  greater  than  unity  as  long  as  carbon  is  present. 
The  assumption  that  the  reaction  goes  completely  to  the  right  results  in 
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an  overestimation  of  the  detonation  rate  since  the  heat  absorption  is 
overcompensated  for  by  an  increase  in  the  number  of  moles  of  gaseous 
products. 

As  a  first  approximation  results  may  be  obtained  using  the 
follovdng  arbitrary  decomposition  equations  for  explosives  of  the 
general  formula,  CqHj.OqN^. 
Case  A:  (q  •►  r/2  ^  S) 

C  HyOgN^  =  r/2  H2O  +  (s  -  q  -  r/2)C02  +  (2q  -  s  +  r/2)C0  +  t/2  Ng 

Case  B:  (q  ••■  r/2  >-   S) 

C  H  OgN^  =  r/2  H2O  +  (s  -  r/2)C0  (q  -  s  +  r/2)C  +  t/2  N2 

These  equations  ejq^ress  the  assumption  that  all  the  hydrogen 
is  converted  to  water,  that  the  remaining  oxygen  then  converts  carbon 
to  carbon  monoxide  and  that  additional  oxygen,  if  any,  converts  carbon 
monoxide  to  carbon  dioxide.  The  results  of  a  calculation  using  the 
arbitrary  product  composition  given  by  the  above  equations  are  shown  in 
Table  III. 

If  one  assumes  that  the  detonation  products  are  in  chemical 
equilibrium  at  the  final  temperature  pressure  and  density,  then  the 
equilibrium  relations  together  with  the  stoichiometric  restrictions 
provide  the  necessary  relations  for  the  calculation  of  the  product 
composition.  The  thermodynamic  equilibrium  constants  for  the  significant 
equilibria  can  be  computed  from  the  function  -  (F°-Eq)/T  for  each  of 
the  components,  together  with  ^Eq  for  the  equilibria  of  interest. 
Tables  of  these  quantities  for  "A^t   CO,  CO2,  H2O,  N2,  OH,  O2,  and  NO 
and  for  the  following  equilibria  over  the  temperature  interval,  2000°  - 
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5000°K.  have  been  published  (30). 

H2  +  CO2  =  CO  +  H2O 

2C0  =  CO2  +  C 

HgO  =  ^H2  ♦  OH 

2H2O  =  2H2  +  O2 

H2O  •»•  ^Ng  =  NO  +  H2 
The  problem  is  conveniently  approached  by  the  consideration  of 
the  work  content  A  of  M  grams  of  a  mixture  of  gases  at  a  large  volume 
V°  and  temperature  T-i .  From  the  thermodynamic  relations,  we  have 
A(V°,  Ti)  =  In^yul   +  RTi  InTrCpi)^^  -  nRTi      (54) 

o 

where  n^   is  the  number  of  moles  of  the  i-th  constituent  and  y**  4  and 

p^  are  its  chemical  potential  and  partial  pressure,  respectively,  in 
the  standard  state  (pure  i-th  component  at  T^  and  unit  pressure). 
n  is  the  total  number  of  moles. 
Equation  34  can  be  written 

A(V<>,  T],)  =  Xni^J  +  RTi  [n  In  p/n  4  Jn^^  In  n^  -  n](5$) 

The  work  contait  at  a  volume  V^  is  obtained  by 

A(Vi,  T^)  =  A(VO,  Ti)   -  M  j^'jpdV  (56) 

To  evaluate  the  integral  requires  the  introduction  of  an  equation  of 
state.  Since  calculations  reported  in  a  later  section  have  been  based 
on  the  Halford-Wilson-Kistiakowsky  equation,  Equation  33*  this  will  be 
used.  To  simplify  the  resulting  expression  it  is  convenient  to  use  the 
fxinction  F(x)  as  defined  in  Equation  33b,  .  Performing  the  indicated 
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Integration  and  simplifying  one  obtains 

A(V;l'^i)  =  Jn^/^l   +  RTi  {n^  In  n^  4  nln  [  P3_/nF(x;i^)] 

n(e^^-.l)  /p-n)  (57) 

The  chemical  potential  is  given  by 

-^  =  (  aA/ ani)j^y^„  (58) 

and  therefore 

^1'  ^X'¥  RTi  [  In  n^  +  ln(Pi/n)  -f  2.303  /I  (x^)  ♦ 

2.303  n  ki  r5(x)/k]  (59) 

where    Q  (x)  «  0.4343  (e  ^^  -  1)/^  -  log  F(x)  (60) 

fz  (^)  =  0.4343  X  e?^  (61) 

These  functions  have  been  tabulated  with  x  as  the  argument  (31)*  The 
chemical  reaction  may  be  represented  by 

ZW^i  x^  =  0  (62) 

where  the  mole  number  W^j  for  the  i-th  constituent  in  the  J-th  reaction 
is  taken  positive  for  products  and  negative  for  reactants.  Thus 
Equation  59  niay  be  written 

^ji>"i  ■   ^ji^i  "*'^^l[^ji  1"  "i  •»  2:Wj^  (In  (Pj^/n)  ^ 
2.303  ri(xi)  +  2.303  n  ZMj^  ki  r2(x)/k]     (63) 

Now,  since  at  equilibrium  Z'Wjj^><^  =  0  and  by  definition  Sw^^^yU^  = 
-  RT  In  Kp  we  may  write 

log  K   =  ZWji  log  n^  -  log  Gj  (64) 

where  Kp  is  the  thermodynamic  equilibrium  constant  for  the  j-th  reaction 


and  the  qiiantity  G.  is  defined  by  the  equation 

-log  Gj  =     [logCPj^/n)  +  Q  (x^)]  ZWj^  +  n  P2  {x^) 

ZWji  k^/k  (65) 

Thus  one  can  now  write  for  the  equilibrium  constant  in  terms  o£ 


moles  (K.) 
J 


'^i  =  "pj  °j '  ■^^°i^ 


"ji 


(66) 


If  there  are  1  constituents  in  the  mixture,  composed  of  m  different 
atomic  species,  there  will  be  (1  -  m)  relations  of  the  form  of 
Equation  66  and  m  stoichiometric  relations.  A  convenient  iterative 
procedure  and  the  tabulated  functions  necessary  for  the  calculation  of 
the  equilibrium  composition  of  the  reaction  products  are  available  in 
reference  (31).  The  equilibria  usually  considered  include  the  five 
equations  given  above  as  well  as 

3H2  f  CO  =  CH^  +  H2O 
and  3/2  H2  +  i  N2  =  NH3 

A  comparison  of  the  final  composition  of  some  typical  ejqplosives 
based  upon  the  arbitrary  decomposition  equation  and  the  equilibrium 
considerations  is  given  in  Table  II  below.  The  composition  is  given 
as  the  number  of  moles  of  each  constituent  per  mole  of  explosive. 
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Table  II 


TNT  P, 

3  =  1.0 

Tetryl  Pq  "  ^'^ 

Cyclonite  ft  =  1.3 

Arbit. 

Equil. 

Arbit. 

Lquil. 

Arbit. 

EquiJ 

«2 

0 

0.09 

0 

0.08 

0 

0.02 

CO2 

0 

0.18 

0 

0.24 

0 

0.16 

CO 

3.5 

3.38 

5.5 

5.27 

3.0 

2.78 

H2O 

2.5 

2.26 

2.5 

2.24 

3.0 

2.87 

N2 

1.5 

1.46 

2.5 

2,45 

3.0 

2.96 

iO 

0 

0 

0 

0 

0 

0.02 

CH4 

0 

0.01 

0 

'0.01 

0 

0 

^«3 

0 

0.09 

0 

0.10 

0 

0.07 

c 

3.5 

3.43 

1.5 

1.48 

0 

0.06 

DH 

0 

0 

0 

0.01 

0 

0.01 

Comparison  of  product  composition,  based  upon  arbitrary  and  equilibrium 
relations. 
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Rli^ULTS  OF  im,   CALCULATIONS 

Table  III  is  an  example  of  the  results  obtained  by  Brinkley 
and  Wilson  (31)  using  the  method  of  approach  discussed  in  the  previous 
section.  The  initial  and  final  densities,  p^  and  p^,  are  in  grams/cubic 
centimeter.  The  detonation  velocity,  D,  and  the  mass  velocity,  u^,  are 
in  meters  per  second.  The  temperature  is  in  degrees  Kelvin  and  the 
pressure  in  kilobars  (lo"  dynes/cm^).  The  composition  is  in  gram  moles 
of  the  individual  constituent  per  gram  mole  of  explosive. 

A  tabulation  of  the  results  obtained  by  Jones  (2?)  using  TNT 
(^Q  -  1.5  is  included  for  comparison.  Jones'  final  composition  gives 
larger  values  for  CO2  and  CH,  at  the  expense  of  H2O  and  CO.  The  avail- 
able experimental  data  on  detonation  product  composition  indicate  much 
smaller  quantities  of  methane,  however,  it  is  entirely  possible  that 
formation  may  have  occurred  during  the  cooling  process.  The  equilibrium 
of  reactions  involving  CO  and  CO2  are  also  sensitive  to  temperature  changes, 
On  the  whole,  it  is  reasonable  to  suppose  that  the  compositions  given 
by  the  Brinkley-Wilson  procedure  are  a  closer  approach  to  the  actual 
conditions.  This  is  also  suggested  by  the  better  agreement  with  the 
observed  (32)  detonation  velocity  of  6620  meter/sec.  for  an  initial 
density  of  1.5, 

The  right-hand  column  contains  results  based  on  the  arbitrary 
decomposition  equations.  It  appears  that  the  detonation  velocity,  which 
is  the  only  property  that  can  be  measured  accurately,  is  compsiratively 
insensitive  to  the  composition.  Although  the  three  methods  give  results 
of  the  same  order  of  magnitude,  it  is  felt  that  the  equilibrium  approach 
of  Brinkley  and  Wilson  is  preferable. 
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A  comparison  of  observed  and  calculated  detonation  velocities 
of  some  common  explosives  is  shown  in  Figures  7,  B  and  9.  The  observed 
velocities  are  those  reported  by  the  l!i3q)losives  Research  Laboratory  (32), 
Bruceton,  Pa.  The  calculated  values  are  based  on  the  method  of  Brinkley 
and  Wilson  vdth  consideration  of  the  seven  primary  equilibria.  The 
overall  agreement  is  considered  to  be  very  good. 

Deviations  may  be  due  to  errors  in  the  observed  quantities, 
however,  these  are  usually  accurate  to  within  3  -  k%,   or  to  errors  in 
the  heat  of  formation  of  the  explosive  or  to  inadequacies  of  the  method 
due  to  lack  of  accurate  data  in  regard  to  covoliome  factors  and  equilibrium 
constants  under  the  extreme  conditions. 

The  heats  of  formation  of  the  explosives  considered  in  Figures 
7,  8  and  9  are  based  on  the  values  reported  by  Schmidt  (33)  in  192A. 
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Table  III 
Calculated  Values  of  the  Detonation  Velocity 
and  of  the  Properties  of  the  Detonation  Wave 

TNT 


« 

it* 

ft 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.5 

1.6 

Pi 

1.347 

1.464 

1.578 

I069O 

1.801 

1.911 

2.019 

1.94 

0 

AB70 

5190 

5540 

5920 

6340 

6820 

7290 

7720 

7210 

"l 

1250 

1290 

1330 

1370 

1410 

1460 

1510 

1750 

Tl 

3260 

3270 

3270 

3240 

3210 

3190 

3170 

3460 

2880 

Pi 

61.0 

73.6 

88.6 

105  oO 

125.4 

151.6 

177.8 

205 

173 

Composition 

H2 

.09 

.05 

.02 

oOl 

0 

0 

0 

0 

0 

CO2 

.18 

.21 

.25 

.31 

.38 

.46 

.54 

2.64 

0 

CO 

3.38 

3.22 

3.08 

2.92 

2.75 

2.58 

2.42 

.06 

3.5 

H2O 

2.26 

2.36 

2.42 

2.46 

2.48 

2.48 

2.48 

1.10 

2.5 

N2 

1.A6 

1.47 

1.48 

1.49 

1.49 

1.48 

1.48 

1.50 

1.5 

NO 

0 

0 

0 

0 

.01 

.02 

.02 

0 

0 

CH4 

.01 

0 

0 

0 

0 

0 

0 

.70 

0 

NH3 

.09 

.06 

.04 

.02 

.01 

.01 

.01 

0 

0 

C 

3.43 

3.57 

3.67 

2.77 

3.87 

3.96 

4.04 

3.84 

3.5 

*  Based  on  Jones'  equation  of  state. 
**  Based  on  arbitrary  decomposition  equations. 
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APPi^DIX  A 

RtXA'flQN  OF  DETONATION  VELOCITY  TO  c  •►  u  AT  VARIOUS  POINTS 
ON  THE  RANKINE-HUGONIOT  CURVE 

In  a  previous  section,  it  was  shown  that  at  the  point  B, 
Figure  3>  the  detonation  velocity  was  equal  to  c  •(■  u  in  accordance  with 
the  Chapman-Jouguet  condition.     It  is  of  interest  to  determine  the 
relation  at  points  other  than  B.     Figure  3  is  reproduced  below  for 
convenient  reference* 


R-H  curve 


Figure  3 
Generalized  Rankine-Hugoniot  Curve 


It  is  well  established  that  the  sonic  velocity  is  related  to 
the  slope  of  an  isentropic  PV  curve  by 

-s = ■?  (-  a 

From  Equations  4.  and  5>  we  have  the  relation 

^   1'     1  V  -V, 
*o  "1 

PqVq  being  evaluated  at  the  point  A  and  PiV^  on  the  R-H  curve.  From 

2 
these  two  equations,  it  appears  that  the  difference  between  o-^   and 


-  53  - 


•3 

(D-Uj^)     at  various  points  on  the  R-H  cvunre  is 


(1) 


The  slopes  of  the  isentrope  and  R-H  curves  at  any  point  (P,  V)  are 
related  by  the  equation 


(3v)g°(l?)/[(as),(|f]pi  a 


where  the  subscript  R  indicates  differentiation  along  the  R-H  curve. 
Rearranging  and  introducing  the  Maxwell  relation f-:^]  =  -  (-Trr) 
we  have 


'  av/s    r^vja  ■      lav 


7p  lav,p 


(2) 


R 


From  the  energy  Equation  6  for  points  on  the  R-H  curve,  we  have 

(1^1  "^(v»i)(^)  -i(W 


(- 


Substituting  for  f   1|  its  equivalent  given  by  the  first  and  second 


laws  of  thermodynamics  yields 


asi\     ^  V!i 


avi/ 


R 


2T. 


Vo-^l 


9P] 


3V 


1/R 


(3) 


The  sign  of  quantity  on  the  left-hand  side  of  this  equation  depends  upon 
the  relative  slopes  of  the  line  joining  the  initial  and  final  states  and 
the  tangent  to  the  curve  at  the  point  representing  the  final  state.  Two 
extremes  are  shown  in  the  figure  below.  Ivote  that  at  point  B  the  two 
lines  coincide. 
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hh 


R-H  curve 


Figure  10 
Slope  Relations  on  the  R-H  Curves. 

If  the  R-H  curve  is  concave  upward!  2r^*  ^^^   difference  in  slopes  is 
negative  above  B  and  positive  below  B.  It  is  certainly  reasonable  to 
e:q)€ct  that  the  R-H  curve  should  have  this  type  cux^ature  for  the 
majority  of  normal  fluids.  Solving  Equation  3  for  the  quantity  within 
the  square  brackets  and  subtracting  from  Equation  2,   evaluated  at  the 
point  (P^V^)  permits  the  velocity  difference  written  in  Equation  1  to 
be  expressed  in  terms  of  variables  of  known  sign.  That  is 


•;— ^'=-(^Ue;.'(^) 


(4) 


The  terms  within  the  square  brackets  are  always  positive  as  is  vf ; 
therefore  the  right-hand  side  is  positive  for  points  on  the  R-H  curve 
above  B  (where  (  93-^/  9V^)rj  is  negative)  and  negative  for  points  below 
B,  that  is  where  (  9Sn/  3Vi)p  is  positive.  Hence,  as  was  to  be  shown, 
D  is  less  than  Cj^fu^^  when  the  final  state  (PiV]_)  is  above  B  and  is  greater 
than  d^fuj^  when  the  final  state  is  below  B, 
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APPENDIX  B 
SECOiPARY  EFFLCTS 

The  theoretical  treatment  of  the  detonation  of  solid  explosives 
was  most  conveniently  accomplished  by  neglecting  certain  effects  such  as 
charge  diameter  and  unstable  detonation  velocities  that  must  be  considered 
in  any  practical  application  of  the  theory.  Fundamental  investigations 
of  these  effects  and  their  ramifications  have  been  made  and  several 
theories  have  been  pzx>posed  in  an  effort  to  explain  the  observed 
phenomena.  Since  a  detailed  description  of  the  proposed  theories  sheds 
little  light  on  the  fundamentals  of  the  detonation  process^  it  is  believed 
advisable,  for  the  purpose  at  hand,  to  regard  them  as  "second-order" 
effects  and  to  devote  to  such  effects  a  very  limited  treatment. 

The  effect  of  cnarge  diameter  on  the  detonation  velocity  has 
been  observed  since  the  earliest  •jq)eriments  were  conducted.  This  effect 
is  observed  as  a  decrease  in  the  detonation  velocity  when  an  explosive  is 
confined  in  a  tube  of  small  diauneter.  Fortunately,  the  variation  ceases 
when  a  certain  reasonable  diameter  is  reached,  thus  the  investigation  of 
other  parameters  can  be  carried  out  without  the  introduction  of  an 
additional  variable.  By  the  same  token,  few  investigators  have  under- 
taicen  detailed  studies  of  its  cause  and  until  recently  a  generally 
accepted  explanation  was  not  available.  A  theory  proposed  by  Jones  (34) 
accounts  for  the  observed  variation  on  the  basis  of  lateral  losses. 
During  the  finite  reaction  time,  the  high  pressure  with  the  reaction  zone 
is  assumed  to  cause  an  appreciable  lateral  e^qiansion  with  a  resultant 
drop  in  pressure,  temperature  and  detonation  velocity.  In  other  words. 
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not  all  of  the  energy  of  the  e^qslosive  is  available  to  maintedn  the  ideal 
detonation  velocity  due  to  lateral  dissipation.  It  can  be  shown  by  a 
detailed  mathematical  analysis  that  the  effect  becomes  significant  when 
the  radius  of  the  explosive  charge  is  of  the  same  order  of  magnitude  as 
the  reaction  zone. 

Another  phenomenon  observed  frequently  in  the  detonation  of 
solid  e:q)losives  is  the  transitioi  from  an  unstable  detonation  velocity 
(one  that  is  either  too  high  or  too  low)  to  the  stable  value.  Whenever 
an  explosive  charge  is  initiated  by  a  priming  charge  having  a  different 
detonation  rate,  this  transition  must  occur.  A  theory  of  the  building-up 
rate  as  a  function  of  the  reaction  zone  length  has  been  proposed  (7)  and 
leads  to  the  conclusion  that  the  ideal  velocity  is  approached  asymtotically, 
half  the  difference  between  the  initial  and  final  value  being  covered 
in  a  few  reaction  zone  lengths. 

If  a  detonation  wave  experiences  such  severe  losses  that  the 
wave  can  no  longer  propagate  at  the  ideal  velocity,  the  phenomenon 
referred  to  as  "fading"  is  observed.  As  in  the  case  of  transition  from 
one  velocitj^  to  another,  the  length  of  the  reaction  zone  is  a  critical 
feature  and  theoretical  treatments  of  the  phenomenon  are  based  upon  the 
cumulative  effects  resulting  from  increased  reaction  zone  lengths. 

In  summary,  it  might  be  observed  that  if  one  possessed  sufficient 
information  concerning  the  chemical  reaction  kinetics  to  predict  the 
length  of  the  reaction  zone  under  any  conditions  of  pressure,  temperature, 
etc.,  a  more  rational  approach  to  many  of  the  technical  problems  con- 
fronting the  ejqslosives  engineer  could  be  made. 
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